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The  final  (minimum  radius)  ^te  variables  of  a cylindrically  imploded  rotating  liquid  liner 
are  computed  as  functions  of  Pf,  the  maximum  compression  achieved,  and  Uaj.,  the  velocity  the 
liner  would  attain  if  allowed  to  expand  without  restraint.  For  each  choice  of  Pf  and  u_,  the 
rotational  speed  is  chosen  to  just  stabilize  the  Rayleigh-Taylor  modes  on  the  inside  surwce.  The 
acceleration  of  the  inner  surface  is  largest  prior  to  turnaround,  so  that  the  rotational  speed  required 
for  stabilization  is  close  to  that  found  for  an  equivalent  incompressible  liner.  Near  turnaround 
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20.  Abstract  (Continued) 

the  inner  portion  of  the  liner  becomes  significantly  compressed,  making  the  efficiency  with  which 
the  payload  plasma  is  heated  considerably  less  than  that  for  an  incompressible  liner.  The  liner 
compression  produced  at  tumaroui^  alters  the  implosion  dynamics  and  creates  a pressure  pulse 
propagating  outward  analogous  to  ^*water  hammer.'**^ An  important  favorable  effect  of  compres- 
sion is  to  extend  the  dwell  period,  during  which  thermonuclear  reaction  rates  are  maximal.  The 
dependence  of  compression  efficiency  and  of  Q,  defined  as  the  ratio  of  thermonuclear  yield  to 
total  system  energy,  on  pf  and  u^o  ate  displayed  for  the  parameter  range  of  experimental  interest. 
These  results  are  then  i|s^d  to  scalje  the  system  size  in  terms  of  allowable  mechanical  stress  to 
indicate  optimal  peak^operating  pi|essures.  The  sidewall  and  endwall  pressure  loads  are  calculated 
for  the  same  parameter  range. 
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ROTATING  IMPLODING  LIQUID  LINERS 


, 

I.  Introduction 

The  Linus  program  at  the  Naval  Research  Laboratory  is  aimed  at 
investigating  the  possibility  of  developing  controlled  thermonuclear 
fusion  by  adiabatically  compressing  plasma  by  controlled  liner 
implosion. The  concept  is  based  on  the  use  of  liquid  metal 

C2*5l 

liners  (rather  than  the  solid  liners  used  in  early  experiments^  ) 

stabilized  by  rotation^^'  imploded  hydraulically  by  high-pressure  I 

gas-driven  pistons.  Working  with  a liner  which  implodes  cyclically. 


rather  than  disintegrating  at  each  shot,  permits  recovery  of  much  of 
the  driving  energy.  This  means  that  the  system  can  operate  at  much 
lower  values  of  Q,  here  defined  as  the  ratio  of  the  fusion  yield  to 
the  total  energy  imparted  to  the  liner.  Since  the  volume  and  the 
total  energy  of  a reactor  system  scale  roughly  as  the  square  of  Q, 
substantial  reduction  in  reactor  size  becomes  possible  in  principle. 

The  adoption  of  a "captive"  (piston-driven)  liquid  liner 
design^^^  eliminates  the  outer  free  surface,  which  if  present  would 
necessarily  exhibit  Rayleigh-Taylor  instability  near  the  outer 
turning  point  (i.e.,  at  the  onset  of  the  implosion  and  near  the  end 
of  the  re-expansion  phase).  Thus,  exploitation  of  the  liquid  pro- 
perties is  essential  to  the  operation  of  a repetitively  pulsed  liner. 

In  this  paper  we  examine  closely  a different  aspect  of  liquid 
liners,  namely,  the  effects  of  the  finite  compressibility  found  in 
real  liquid  metals.  Measured  by  the  simplest  criterion,  the  flow 
Mach  numbers  M (relative  to  c,  the  speed  of  sound  in  the  uncompressed 
material),  these  effects  might  be  expected  to  be  small  or  even 
negligible.  In  liquid  alkali  metals,  for  example,  c **  1.5-7.5  km/sec, 
whereas  the  liner  speeds  attainable  with  hydraulic  compression  tend 
to  be  on  the  order  of  a few  hundred  meters  per  second,  so  that 
M ~ .05-. 2.  (These  values  are  typical  during  most  of  the  Implosion; 
Just  prior  to  turnaround,  speeds  approaching  1 km/sec  are  achieved 


Note:  Manuscript  submitted  January  13,  1978. 
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due  to  convergence  effects.)  But  even  though  the  gross  dynamics  of 
an  implosion  is  practically  incompressible,  the  parameter  of  most 
concern  to  us  - the  amount  of  plasma  compression  - is  greatly 
affected.  For  thermonuclear  applications  the  liner  has  to  be  very 
thick  relative  to  the  inner  radius  at  turnaround;  volumetric 
compressions  ~ 10^,  i.e.,  aspect  ratios  ~ 50,  are  typical.  For  a 
given  value  of  the  outer  radius,  an  average  compression  throughout 
the  liner  of  as  little  as  0.1^  can  double  this  inner  radius  and  reduce 
the  temperature  achieved  by  a factor  ~ 2 2.5 

Another  deleterious  effect  of  compressibility,  known  as  "water 
hammer,"  arises  from  the  nature  of  sound  wave  propagation  in  liquids. 
The  liner  is  stopped  so  abruptly  at  turnaround  by  the  buildup  of 
pressure  in  the  "payload"  (the  mixture  of  plasma  and  magnetic  flux 
being  compressed)  that  a compression  wave  is  launched  outward,  follo- 
wed by  a rarefaction  wave.  This  pressure  pulse  propagates  to  the 
side-  and  endwalls  of  the  system,  producing  stresses  equal  to  a 
significant  fraction  of  the  peak  payload  pressure  and  therefore 
substantially  in  excess  of  the  initial  driving  pressure.  Such  loads 
can  lead  to  failure  of  rigid  structural  members  of  the  device. 

While  carrying  out  numerical  studies  of  these  phenomena,  how- 
ever, we  discovered  that  the  news  is  not  all  bad.  During  and 
immediately  after  turnaround  the  outer  portion  of  a compressible 
liner  continues  to  Implode  and  compress  the  innermost  portion. 

While  this  continues  (roughly  twice  the  sound  transit  time,  defined 
as  the  liner  thickness  d divided  by  c)  the  liner  inner  radius 
expands  at  a slower  speed  than  that  with  which  it  imploded.  This 
prolongs  the  "dwell"  of  the  liner,  thereby  enhancing  Q.  The  effect 
is  naturally  more  pronounced  at  high  Mach  numbers  and  with  thicker 
liners. 

There  is  considerable  uncertainty  about  the  appropriate  form 
for  the  equation  of  state  of  liquid  metals  at  high  pressures  and 
temperatures.  Earlier  treatments  concentrated  on  the  behavior  of 
solid  liners.  Knoepfel^  summarized  the  available  measurements  and 

[91 

theoretical  models  and  quoted  a quasi-emplrical  formula  of  Somon 
for  pressures  less  than  1 Mbar.  We  use  this  form  below  in  our 
calculations  even  though  it  was  derived  for  solid  state.  For 
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hydrodynamic  applications  the  exact  functional  form  turns  out  not  to  be 
important,  since  at  Linus  pressures  “ 10^  bar)  the 


velocity  of  sound  deviates  only  slightly  from  its  zero-pressure 
value. 

The  model  studied  in  the  present  paper  differs  in  significant 
respects  from  those  treated  in  previous  work.  Earlier  treatments 
all  left  rotation  out  of  the  picture;  in  Linus,  by  contrast,  when 
Rayleigh -Taylor  instability  is  fully  suppressed,  roughly  half  of 
the  system  energy  ends  up  as  rotational  energy.  In  addition, 
previous  treatments  are  inappropriate  for  Linus -type  devices  in 
other  ways.  For  example,  Lehner  ^ al^^*^^  studied  magnetic  flux 
compression  produced  by  a slab  liner,  which  relates  closely  to  the 
cylindrical  case  only  if  the  liner  thickness  d is  much  less  than 
the  final  radius  r^.  They  introduced  the  concept  of  the  effective 
thickness  d*,  defined  as  the  thickness  an  incompressible  liner 
with  the  same  density  o would  have  to  have  to  produce  the  same  final 
payload  pressure  p^  for  a given  initial  liner  displacement  r^  » r^ 
and  implosion  speed  u^^  * H^c.  Only  that  portion  of  the  liner  within 
a distance  d*  from  the  inner  surface  can  participate  effectively 
in  the  compression. 

The  resfults  they  obtained  may  be  summarized  as  follows.  The 
effective  thickness  is  given  by 


ti  ^Pcu^^ 


p i (1-KXM  )i 
I \2  1 


(1) 


M. 


where  for  magnetic  flux  compression  (Y  * 2) 


"i  ■ 


(2) 


is  the  initial  pressure  of  the  payload  and  a 1.5  is  a phenomenolo- 
gical constant.  For  thick  liners  the  final  pressure  achieved  is 
given  by 




Pf/PC‘ 


M^(l-KlM^). 


(3) 
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By  replacing  Eq.  (2)  with  the  two-dimensional  form, 


Pi  = Pf  (r^/r.)'' 


we  can  define  the  compression  efficiency  for  an  equivalent  cylindri- 
cal geometry  according  to 

fl)  ""  Pf 

* “ (2TT  r^d^)  (tP  ^5) 

where  ratio  of  the  Initial  cross-sectional 

area  of  the  liner  to  that  of  the  final  payload  volume.  Substitution 
of  (3)  in  (5)  yields 


2(1  +aM^) 


bf-^Mi 


The  efficiency  cannot  exceed  unity,  of  course;  as  0,  Eq.  (6) 

must  go  over  to  g = 1 (the  incompressible  result).  In  the  high 
Mach  number  case,  i.e.,  when  Eq.  (3)  is 

replaced  by 


Cj  - M^(l4aM^)  (d^/d*)2  « \ M^b^2 


whence  from  (5) 


(2M^)- 


Somon  and  Jablon  extended  these  results  to  the  case  of 

supersonic  (M^»l)  cylindrical  liners.  Their  numerical  simulations 
showed  that  the  maximum  pressure  is  given  by 


- 0.89  5 
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(9) 


and  hence  the  compression  efficiency  is 


(10) 


Other  quantities  of  interest  can  be  found  from  and  Eq. 

Robson  applied  the  concept  of  the  effective  liner  thickness 
to  cylindrical  geometry  for  the  low  Mach  number  case  and  final 
pressures^^  ^ 0.05.  The  results  are  rather  Insensitive  to  the  choice 
of  the  "effective  outer  radius"  r*.  They  show  that  the  ratio  of 
p^  to  p^,  the  final  pressure  the  liner  would  reach  if  it  were  incom- 
pressible, is  a decreasing  function  of  both  and  b^ (parameters 

defined  analogously  to  Cc  and  b using  p Instead  of  p,  and  r , the 

r r o to 

incompressible  final  state  radius,  instead  of  Our  results  in 

the  present  paper,  which  are  restricted  to  ^ 0.05  for  numerical 
reasons  discussed  below,  appear  to  be  roughly  consistent  with 
Robson's. 

Schaffer,  using  the  same  equation  of  state  and  the  same  fluid 
equations  of  motion  (except  for  the  absence  of  rotation)  as  ours, 
performed  a series  of  numerical  experiments  using  an  implicit 
finite-difference  algorithm.  He  found^^^^  that  a good  fit  to  the 
results  could  be  expressed  in  the  form 

e Cl  + Cf(b^/F)2]’^ 

where  ^ 

'jl+0.65  ^n(b^/20)  b^  >20 

Over  a considerable  range  of  the  parameters  the  effective  radius  r 
satisfied  26. 

In  order  to  understand  better  the  effects  of  compressibility  and 

to  develop  scaling  laws  of  a nature  sufficiently  general  to  aid  in 

[i41 

designing  and  interpreting  liquid  liner  experiments  > we  have 
adopted  the  following  simple  model  (cf.  Fig.  1).  A cylindrical  liner 
is  described  using  the  Euler  equation  for  the  radial  aximuthal 
velocity  components,  together  with  the  equation  of  state  given  by 


fA.8 


b^  20 


(11) 


(12) 


Somon 


[9] 


Both  liner  surfaces  are  assumed  to  be  free.  The  pressure 


5 


at  the  inner  surface  is  derived  from  the  adiabatic  law  for  the  payload 


Pl  ■ "o 

where  Y is  the  adiabatic  index  (usually  taken  to  be  5/3)  and  is 
the  pressure  at  the  inner  radius  r^.  The  outer  pressure  is  taken  to 
be  constant  (specifically,  zero  for  the  calculations  reported  here). 

Thus  the  liner  and  plasma  are  taken  to  behave  as  ideal  fluids  no 
energy  dissipation  or  other  loss  mechanisms  are  modeled). 

By  neglecting  the  driving  pressure  and  the  boundary  conditions 
appropriate  to  realistic  driving  pistons,  we  greatly  reduce  the 
complexity  of  the  model,  without  lessening  its  applicability  to 
the  part  of  the  trajectory  in  the  immediate  vicinity  of  turnaround. 
As  in  all  liner  implosion  problems,  the  most  interesting  effects 
which  we  wish  to  concentrate  on  are  observed  only  near  turnaround 
anyway.  The  liner  rotation  is  then  chosen  in  such  a fashion  that 
the  criterion  for  rotational  stabilization  of  Rayleigh-Taylor  modes 
on  the  inner  surface  is  marginally  satisfied  during  the  implosion: 


where  Vj^  is  the  azimuthal  velocity  at  the  inner  surface,  [it  turns 
out  that  the  maximum  in  g(t)  always  occurs  just  prior  to  turnaround, 
where  equality  therefore  holds  in  Eq.  (1^).] 

Using  this  model,  we  initialize  the  liner  trajectory  at  some 
point  shortly  before  turnaround  and  solve  the  equations  of  motion  in 
finite-dlfferenced  form  to  find  the  behavior  of  the  liner  as  a 
function  of  time.  The  fractions  of  total  energy  at  turnaround  in 
the  form  of  payload  compression,  rotation,  liner  internal  energy  and 
radial  kinetic  energy,  the  value  of  Q,  and  the  amplitude  of  the 
pressure  pulse  at  the  outer  radius  are  found  for  various  initial 
conditions.  The  results  are  summarized  in  the  form  of  plots  versus 
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and  the  limiting  form  of  u^/c  at  asymptotically  large  radius, 
then  crossplotted  to  reveal  their  dependence  on  effective  driving 
pressure  and  compression  ratio. 

The  plan  of  this  paper  is  as  follows.  In  Section  II  we  review 
the  dynamics  of  incompressible  liners,  then  recast  the  problem  in 
dimensionless  variables.  The  numerical  techniques  used  are  described 
in  considerable  detail,  especially  the  way  in  which  knowledge  of  the 
incompressible  dynamics  is  used  to  initialize  the  calculations.  The 
results  of  the  calculations  are  then  presented  in  Section  III,  with 
particular  attention  to  the  wave  effects  (water  hammer  and  dwell 
enhancement)  mentioned  above.  We  summarize  our  conclusions  in 
Section  IV. 

II.  Formulation  of  the  Problem 
A.  Review  of  Incompressible  Liner  Dynamics 

We  begin  our  treatment  of  compressibility  effects  in  imploding 
cylindrical  liquid  liners  by  reviewing  the  most  important  features 
of  the  motion  of  incompressible  liners.  Although  the  latter 
constitute  an  idealization  quite  unattainable  in  operating  regimes 
of  reactor  interest,  their  study  is  valuable  for  two  reasons.  First, 
they  are  easier  to  analyze,  and  for  some  purposes  the  effects  of 
compressibility  can  be  ignored  (but  it  is  important  to  know  when  they 
cannot).  Second,  at  sufficiently  large  radii  any  freely  imploding 
or  expanding  liner  behaves  as  if  it  is  incompressible.  Therefore, 
the  entire  dynamic  trajectory  can  be  specified  simply  by  writing 
dovm  the  parameters  which  specify  its  asymptotic  incompressible  form. 

As  Somon^^^^  has  shown,  given  the  equations  of  state  of  liner 
and  payload,  two  dimensionless  parameters  are  required  for  this  pur- 
pose in  the  absence  of  rotation.  Inclusion  of  an  arbitrary  free 
vortex  rotation  (azimuthal  velocity  inversely  proportional  to  radius) 
increases  the  number  to  three.  That  is,  any  physical  quantity  q 
associated  with  a particular  liner  trajectory  is  determined  as  a 
function  of  time  according  to  a law  of  the  form 


q*qCt;p,p^p;S,S,S], 

12  3 12  3 


where  the  p^  are  dimensionless  parameters  labeling  the  trajectory, 
and  the  are  dimensionless  quantities  which  set  the  scales  of 
length,  time  and  mass. 

For  an  incompressible  liner  it  is  natural  to  choose  as  scale 
factors  the  liner  density  P,  the  minimum  inside  radius  r^,  and  the 
radial  "velocity  at  infinity"  u^,.  The  latter  is  the  uniform  coasting 
speed  approached  in  the  absence  of  external  pressure  at  radii  so 
large  that  the  liner  is  thin  and  the  back  pressure  of  the  compressed 
payload  and  the  dynamic  pressure  within  the  liner  resulting  from 
cylindrical  convergence  are  negligible.  The  dimensionless  parameters 
are  conveniently  chosen  as  b,  the  ratio  of  liner  thickness  to  inner 
radius  at  turnaround;  the  angular  velocity  of  the  inner  surface 

at  turnaround,  scaled  by  u^/r^;  and  P^,  the  payload  pressure  p^  at 
turnaround,  divided  by  pu^^.  The  zero  of  time  is  either  taken  as  the 
instant  of  turnaround,  or  as  the  time  when  we  begin  following  the 
trajectory,  in  which  case  turnaround  occurs  at  t = t^.  (Throughout 
this  paper,  the  zero  subscript  denotes  quantities  associated  with 
incompressible  liner  motion,  evaluated  at  turnaround.  The  corres- 
ponding values  for  compressible  liners  are  labeled  with  the  sub- 
script f,  for  "final.")  This  choice  of  dimensionless  parameters  and 
scale  factors  has  the  advantage  of  relating  the  dynamical  trajectory 
to  the  properties  of  the  system  at  turnaround,  i.e.,  those  which 
determine  the  thermonuclear  properties  of  the  system.  It  may,  how- 
ever, not  be  optimum  for  purposes  of  experimental  design. 

The  equation  of  motion  is  the  radial  component  of  the  force  law. 


+ ^ I?  - 0, 


where  u and  v are  the  radial  and  azimuthal  components  of  velocity 
and  p is  the  pressure.  Conservation  of  angular  momentum  is  expressed 
by  the  azimuthal  component  of  the  force  law. 
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5v  , fiv  . uv 

TT  + u jT"  — 
6c  ar  r 


0. 


(17) 


The  equation  of  state  is  the  condition  of  incompressibility, 

(ru)  = 0 (18) 

By  virtue  of  (l8),  the  radial  velocity  field  in  the  liner  is 
inversely  proportional  to  r, 

u(r,t)  = r (t)  u (t)/r  (19) 

1 1 

(subscripts  1 and  2 label  quantities  associated  with  the  inner  and 
outer  surfaces  of  the  liner,  respectively),  and  the  area  between  the 
inner  radius  and  a radial  displacement  comoving  with  a liner  fluid 
element  is  a constant, 

r^  - r 2 = r^(0)  - r ®(0)..  (20) 

1 1 

Solution  of  (17)  yields  the  result  that  v is  a linear  combination 
of  a free  vortex  and  a rigid  body  rotation.  We  omit  the  latter  for 
simplicity  (this  is  a good  approximation  if  the  liner,  which  starts 
from  a state  of  rigid  rotation,  is  thin  at  that  time).  The  azimuthal 
velocity  field  of  the  liner  can  now  be  written 

v(r,t)  = n^u^  r^/r.  (21) 

Substitution  of  (19)  and  (21)  in  (I6)  enables  us  to  calculate 
the  pressure  profile 

p(r,t)  » p + [p  -p  -6p  (1-r  ^/r  ^)]  ^n{r^/r 
121  12  1 

/ in(r  2/r  2)  + 6p  (1-r  2/r2),  (22) 

21  1 
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where 


P (u^  + v^), 
1 1 


The  pressure  p on  the  outer  surface  is  assumed  to  vanish,  so  that 
2 

we  do  not  have  to  model  the  mechanism  that  produces  the  Implosion. 
With  this  external  pressure  Is  In 

fact  negligible  for  all  states  In  the  vicinity  of  turnaround. 
Equation  (22)  has  a maximum  at  r * r^,  defined  by 


r ^ 6p  £n(r  ^/r 

_SL.  „ 2 ^ (2U) 

r ^ p + 6p  (1-r  2/r 
11  12 

When  6p  Is  sufficiently  large,  r^  ^m  ^2^  there  Is  a maximum 
In  the  pressure  across  the  liner;  then  the  slope  of  the  pressure  at 
r » r^  Is  positive.  From  Eq.  (l6),  this  means  that  the  net  accel- 
eration (Inertial  minus  centrifugal)  on  the  Inner  surface. 


Is  negative.  This  Is  the  criterion  for  rotational  stabilization  of 

c?] 

Raylelgh-Taylor  modes  on  the  Inner  surface  . If  at  some  time  the 

slope  of  p(r)  turns  negative  at  r > r^,  the  motion  becomes  unstable. 

To  estimate  how  large  the  maximum  pressure  p(r  ) Is,  we  can  sub- 

m 

stltute  Eq.  (2h)  In  Eq.  (22)  with  p^  * 0.  The  resulting  expression 
simplifies  If  we  assume  (r^/rg)^«  (r^/rg)^«l,  whereupon 

p(r  ) » p + 5p  * p . (26) 

tn  ^ tnax 

Evidently,  we  can  have  pressure  peaks  In  the  Interior  of  the  liner 
which  exceed  p^  by  a considerable  amount. 

We  assume  that  the  energy  is  conserved  (no  losses  take  place 
through  magnetic  diffusion,  streaming  out  the  ends,  etc.).  The 
pressure  of  the  compressed  payload,  and  therefore  p^,  satisfies  an 


adiabatic  law 

P3^(t)  = P^  Pu„2  (r^yr^)2V,  (27) 

with  ratio  of  specific'  heats  Y satisfying  5/3  ^ Y ^ 2.  This 
corresponds  to  a payload  varying  between  pure  ideal  gas  or 
plasma  (Y  = 5/3)  and  pureaxial  magnetic  flux  (y  “ 2).  By  virtue  of 
Eqs.  (19-22)  and  (26)  it  is  possible  to  specify  the  implosion 
dynamics  totally  once  r^(t)  is  known.  An  equation  for  r^  can  now 
be  written  in  the  form 

(RR  + R^)  -^nd  + b^/R^)  » (r2r2  + Cl  2)  b^/ (R^  + b^) 

o 

+ 2 (28) 

where  R » r /r  and  the  dots  denote  differentiation  with  respect  to 
1 o 

the  dimensionless  time  t * integral  of  Eq.  (28)  is 

/R^R^  + n -)  in(l+b2/R=2)  + 2P  r2(i-y)  /(y-l)-  b^.  (29) 

o o 

This  is  the  dimensionless  form  of  the  statement  of  the  conservation 
of  energy. 


w,  . + w ^ + w » 

kin  rot  pay  tot^ 


where 


w..  _ - i tt  p u 2 r 2 


is  the  total  energy  per  unit  length,  and 


^ ^ P (l+b2/R2), 


* 4-  TT  p u„2  in  (l+b^/R^), 


W - P r2(i-y) 

pay  Y-l  o 


1 


are  the  radial,  rotational  and  payload  energies,  respectively,  per 
unit  length. 

At  turnaround  (tst  =ut/r),R*l  and  R = 0,  and  we  have 
from  Eq.  (29) 


2P. 


(i+b2)  + 


° b2. 


(5^) 


The  compression  efficiency  is 


2P. 


€ = W /W 


(Y-l)b2 


(35) 


We  can  specify  by  demanding  that  at  turnaround  the  system  be 

marginally  stabilized  against  Rayleigh-Taylor  modes  by  rotation. 

The  condition  for  this  to  occur  is  that  Eq.  (1^4-)  vanish  at 

t ■ t , i.e. , 
o'  ’ 


« ('o’  ■"»' 


(3S) 


It  is  easy  to  see  that  if  Eq.  (56)  holds  at  turnaround,  g < 0 for 
all  r^  > r^.  Substituting  Eq.  (36)  in  Eq.  (28)  with  R * 1 yields 

n 2 [in(l+b2)  - b^/d+b^)]  = 2P  . (37). 

o o 


Now  Eq.  (35)  becomes 


Cl  + (V-1)  ^n  U+b^),^  1 

^ ^n  (1+b^)  - b^/d+b^)-* 


-1 


(38) 


Note  that  also  has  dropped  out  of  this  expression  for  e;  this 
result,  that  the  compression  efficiency  depends  only  on  b,  is 

r i6i 

characteristic  of  incompressible  liner  dynamics  . 

In  the  thin  liner  limit  b "*  0,  the  condition  for  rotational 
stabilization  becomes  harder  to  satisfy,  and  e vanishes  as  b^. 

The  efficiency  increases  monotonically  with  b,  and  in  the  thick 
liner  limit  (b  “*  ® ),  e approaches  Y Practically  speaking,  this 
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limit  Is  attained  when  b > 20.  The  fraction  of  the  total  energy 
tied  up  in  rotation  in  a stabilized  liner,  5 “ 1 - Y Is  even 


for  thick  liners. 


B.  Compressible  Dynamics 

The  most  obvious  qualitative  effect  of  compressibility  is  a 
reduction  of  the  payload  compression  achieved  at  the  time  of  turn- 
around, defined  by  dr^/dt  ■ 0.  For  a given  compressible  liner 
trajectory,  asymptotically  identical  at  large  radii  to  that  of  an 
incompressible  liner  of  the  same  dimensions,  radial  and  angular 
velocity  and  density,  the  minimum  radius  actually  reached  will  be 
larger,  r^  > r^.  Some  of  the  energy  that  would  have  gone  into 
payload  compression  goes  instead  into  compressing  liner  material. 

At  the  same  time,  some  energy  will  be  tied  up  in  the  radial  motion 
of  the  outer  portion  of  the  liner,  since  a finite  interval  comparable 
with  the  transit  time  of  sound  is  required  to  communicate  the  infor- 
mation that  the  inner  surface  has  stopped.  In  the  energy  conserva- 
tion equation,  which  now  takes  the  form 


W,  . + W ^ + W + w.  ^ 

kin  rot  pay  int  tot' 


where  is  the  energy  which  goes  into  compressing  liner  material, 

all  the  terms  on  the  left  hand  side  are  in  general  non-vanishing 


at  turnaround. 


An  additional  class  of  phenomena  is  associated  with  compres- 
sional  waves  propagating  in  the  liner.  The  most  dramatic  is  the 

"water  haimner,"  an  abrupt  pressure  peak  which  is  followed  by  a 

rarefaction,  both  expanding  radially  outward  following  turnaround. 
Depending  on  the  external  boundary  conditions,  it  can  reflect  back 
in  and  execute  several  transits  during  the  expansion  phase  of  the 
liner.  Closely  related  to  the  water  hammer  is  an  effect  we  term 
"dwell  enhancement".  For  an  incompressible  liner,  is  symmetric 

about  t “ t^,  and  the  dwell  time  t^,  arbitrarily  defined  as  the 

length  of  time  during  which  ^ exactly  bisected  by 
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turnaround.  When  compressibility  is  accounted  for,  the  liner  is 
seen  to  recede  from  the  turning  point  more  slowly  then  it 
approached.  It  is  in  effect  tamped  by  the  still- imploding  outer 
portion,  and  cannot  expand  freely  until  the  rarefaction  wave 
reflected  from  the  outer  surface  reaches  the  inner  surface. 

To  model  compressibility  in  liquid  metals  realistically,  we 
use  instead  of  Eq.  (l8)  an  equation  of  state  previously  employed  by 

[15] 

Somon" 

P - P (#  - 1).  (UO) 


Properly  speaking,  p should  have  an  additional  term  which  depends  on 

entropy  through  the  temperature  T.  By  adopting  the  isentropic  foim 

p and  0 constant^,  we  relinquish  any  prospect  of  correctly  describing 

shocks,  in  return  for  a considerable  simplification  of  the  numerical 

analysis.  This  point  will  be  further  discussed  below.) 

Equation  i^)  implies  that  the  compressional  energy  density  is 

quadratic  in  P.  The  corresponding  bulk  modulus  of  compression 
dp 

B ■ 0 TT  c*n  be  written 
dc 


B(p)  - B^  + B^'  p. 


where  B^  ' ■ 5 •nd  B^  ■ 0^  is  weakly  dependent  on  p in  the  regime  of 
moderate  compression.  That  is,  Eq.  (Uo)  is  approximately  equivalent 
to  a Griineisen  equation  with  coefficient  B^  ' * 5. 

By  thus  assuming  that  all  fluid  elements  have  the  same  entropy 
and  the  same  specific  angular  momentum  A ^ rv  at  all  times,  we  can 
combine  Eqs.  (I6),  (19)  and  (^4-0)  in  the  form 


5u.  5u  A*^,p,,p  _N 

" a;  - f <5  7 - 


I ^ 


We  also  need  the  continuity  equation. 


It  r (rpu)  - 0. 


Equations (U2) and (45)  can  be  used  to  advance  u and  p,  with  all  the 
other  fluid  quantities  derived  from  these.  The  boundary  conditions 
are  as  before,  Eq.  (2?)  and  = 0. 


It  is  useful  to  write  Eqs.  (42)  and  (45)  in  dimensionless  form, 
using  for  this  purpose  the  uncompressed  density  P,  the  speed  of 
sound  in  the  liner  at  zero  pressure,  c,  defined  by 

= p/P^  (44) 

and  some  approximation  for  the  final  radius  r^.  Since  r^  is  not 
known  a priori,  but  only  as  a result  of  propagating  Eqs.  (42)  and 
(45),  we  begin  by  scaling  with  r^  instead.  Both  are  unambiguously 
determined  by  specifying  the  dynamic  trajectory,  but  r^  has  no 
simple  relation  to  r^.  Thus  we  write 


u'  » u/c;  (45) 

r*  - r/r^;  (46) 

t'  » ct/r^;  (47) 

p'  = p/p;  (48) 

and 

“o  “ A/(r^c).  (49) 


Now  we  have  the  equation 

•^fr-  + ^ Itt  (r'p'u')  = 0;  (50) 

(p'u')+ ^l^r'p'u'2)  + 

(51) 

+ -^rr-  li  + (|  p'  - 2)  p']  - 0, 

where  Eq.  (51)  has  been  rewritten  in  conservative  form.  These  are 
the  basic  equations  of  our  model,  the  numerical  solution  of  which 
provides  the  basis  for  the  present  paper.  Henceforth,  it  is  under- 


1 

1 
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stood  that  we  are  working  with  the  reduced  variables  defined  in 
Eqs.  and  the  primes  will  be  omitted: 


If  + 7 I;  (rpu)  = 0; 


1^  (pu)  - 7 I;:  (rpu^)  + 1;^  p.2)D] 


r 5r 


6r 


(52) 


0.  (53) 


We  now  imagine  that  a liner  has  imploded  to  some  radius  r^^^r^^, 
where  at  time  t»0  we  begin  following  it  by  means  of  Eqs.  (52-53). 

At  some  t .rlier  time,  for  radii  r^^  » r^^,  the  liner  must  certainly 
have  behaved  incompressibly,  and  we  can  characterize  its  trajectory 
by  specifying  b,  and  as  before.  Anticipating  the  Importance 
of  compressibility  effects,  however,  we  employ  instead  the  three 
parameters  M,  and  where 


and 


M - u /c. 


U)  = M.  0 , 
o o' 

o o 


(5^) 

(55) 

(56) 


The  incompressible  thickness  b can  be  derived  from  the  new 
trajectory  parameters  by  Eq.  (3^;),  which  now  becomes 


U)  2 (l+b^)  + 2Q  /^Y-1)  » b^  if. 

o o 


(57) 


The  reason  for  switching  from  b as  a trajectory  label  is  that  at 
t ■ 0,  when  we  initialize,  the  liner  must  already  have  become  some- 
what compressed  [remember  that,  as . shown  by  Eq.  (26),  the  pressures 
within  the  liner  may  become  large  long  before  Pj^  does].  Ignoring 
this  and  taking  r^®  - r^^  to  be  the  same  as  its  asymptotic  value 
would  increase  the  mass  of  the  liner  [according  to  Eq.  (Uo)]  by  an 
amount  dependent  on  the  choice  of  r^.  Under  the  circumstances  b, 
defined  as  the  ratio  of  two  quantities  which  no  longer  have  a simple 
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interpretation  in  terms  of  observables^  becomes  dispensable. 


The  parameters  and  are  still  open  to  the  objection  that 
I their  definition  involves  quantities  associated  with  the  final  radius 

reached  by  an  incompressible  liner.  Therefore,  once  we  have  found 
the  minimum  radius  r^  actually  reached,  we  make  a final  replacement, 
introducing  in  their  stead  the  parameters  and  defined  by 


ID-  * v-/c  = tti  /a, 
t z o t 

(58) 

and 

-2Y 

Cf  “ Pf/P  c2  = Cq  Sf  , 

(59) 

where 

a,  * r,/r 
f f o 

(60) 

Clearly,  uj^  and  Sj,  although  defined  only  implicitly  through 
numerical  integration  of  (52)  and  (53),  are  characteristic  of  the 
trajectory  and,  together  with  M,  specify  it  completely. 

C.  Numerical  Treatment 

At  t * 0,  we  expect  the  functions  u(r),  v(r)  and  p(r)  to  be 
almost  the  same  as  they  would  be  for  an  incompressible  liner, 
[otherwise,  we  have  chosen  a starting  point  too  close  to  turnaround.] 
Yet,  as  we  have  seen,  this  form  of  p(r)  is  inconsistent  with  a 
constant  density  p when  we  adopt  the  equation  of  state  (4o).  For 
reasons  to  be  discussed  shortly,  it  is  impractical  to  choose  r^^  so 
large  that  the  discrepancy  is  negligible.  Instead,  we  proceed  as 
follows.  Write 


u(r,o)  - + u^^\r). 

(61) 

where 

(o)/  ^ • '/ 

u (r)  - r^ri/r 

(62) 

and  we  have  assumed  that 

lu^^^r)!  « 1 u^°^r)l. 

(65) 

17 


Here  ri  is  a constant  to  be  determined.  Substituting  (6l)  in 
Eq.  (53)  and  dropping  first  order  terms  yields 

rir'i  + + m ^ b r 5 

; pi ^ ^ Cp(|  p-2)]  = 0,  (6k) 


where  ri  is  a second  undefined  constant.  Integration  results  in 
quadratic  for  p(r):  • 

P(J  P-2)  -t  [(rir'i  + rl  )=  in  (r^/n^) 


+ (ri^fi  + !U^^)(ri"^  - r"^)]  = -i  [l 
+ (rir'i  + r/^)  ■^n  /^i  ^)  + (fi^  r'l 

+ 00^2)  (n"^  - )],  (65) 

where  r2  is  a third  undefined  constant.  [Note  that  p = 1 at 
r * ra  in  (65)  as  a result  of  the  choice  of  the  integration  con- 
stant.] 

The  three  unknowns  r^,  r^  and  r2  are  the  velocity  and  radial 
acceleration  at  the  inner  surface,  and  the  radius  of  the  outer  sur- 
face, respectively.  They  are  chosen  by  requiring  that  the  pressure 
at  r * ri  have  the  correct  value 

Cl  * 

and  by  requiring  that  the  total  liner  mass  (and  therefore  the  angular 
momentum)  and  the  total  system  energy  have  the  same  values  as  "at 
infinity;" 

r's 

2 ^ r dr  P(r)  - b^,  (67) 

J 

ri 
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I 


rdr  Pvr)[(r.‘-r 


1 -1  (P-l)^l 

2C  r 2(i-Y) 

_£.J: = b^M^. 

Y-1 


(68) 


Equations  (66-68)  were  solved  by  a version  of  Newton's  method,  using 
the  incompressible  values  of  r^,  r^  and  r^  as  the  starting  point. 

The  first  order  correction  u(^)  can  be  obtained  from  Eq.  (52): 


(69) 


r Br  ^ ■ '^bt  ^ “ Br' 

where  the  time  derivative  on  the  right  hand  side  has  to  be  obtained 
by  differentiating  the  implicit  time  dependence  of  r^^.  Similarly, 
first  order  corrections  to  p can  be  obtained  from  Eq.  (53)*  The 
procedure  is  very  involved,  and  in  any  case  cannot  be  expected  to 
yield  a convergent  expansion.  Fortunately,  there  is  no  need  to  go 
beyond  lowest  order. 

In  operation  the  initialization  was  completed  by  specifying  M, 
and  a starting  value  of  the  inner  radius  r^.  The  solutions 
were  propagated  in  time  and  the  function  g(t)  was  monitored  until  u^ 
became  positive,  then  was  incremented  by  some  preassigned  value 


and  the  whole  process  repeated 


After  the  minimum  of  g(t) 
(including  only  times  after  the  grid  ringing  had  ceased)  became  nega- 
tive, the  loop  on  was  restarted  with  a smaller  increment  6C^. 

When  three-place  accuracy  had  been  attained  in  the  crossing  value, 
the  system  was  considered  to  be  marginally  stable,  and  a final  run 
was  carried  out  through  the  expansion  phase  (out  to  some  maximum 
radius  r ~ 5"10  r.)  for  the  purpose  of  calculating  water  hammer 

pressure  histories  and  the  total  thermonuclear  yield.  The  latter  was 

[17] 

calculated  according  to  the  following  formula  for  a 50-50  D-T 
mixture  (njj*n^*0.5n)  as 

Pdt  - '^•3  n- 


t"^  exp  (-19.9^  T"^)erg/cm3 


(70) 


where  T is  in  keV,  and  both  n and  T are  supposed  to  vary  adiabati- 
cally.  Each  calculation  integrated  Eq.  (70)  for  50  different 
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starting  values  of  T,  and  Q was  determined  using  the  one  which  gave 
the  biggest  yield.  Note  that  the  numerical  coefficient  in  (70)  is 
larger  than  that  given  in  ref.(17)  by  (5.0/5.683 (22. 4/17.6)=!. 75.  The 
first  factor  arises  from  optimizing  the  theoretical  fit  to  the 
experimental  cross-section  in  the  energy  range  T ^ 25  keV,  and  the 
second  results  from  taking  account  of  the  energy  produced  by 
tritium  breeding  in  the  lithium  component  of  the  liner. 

As  we  are  interested  in  the  compressional  waves  which  can 

propagate  through  the  liner^  our  code^  called  WAVER,  used  a time- 

centered  explicit  finite  difference  technique  to  solve  Eqs.  (52-55) 

numerically.  An  implicit  technique  would  correctly  describe  the 

gross  features  of  the  implosion,  at  the  same  time  relaxing  the 

r 1 qT 

Courant-Friedrichs-Lewy  (CFL)  timestep  condition  so  that  the 

code  could  run  up  to  two  orders  of  magnitude  faster,  but  it  would  not 

be  accurate  in  describing  wave  propagation  unless  the  time  steps  were 

shortened  and  made  comparable  to  the  CFL  limit.  Because  of  the 

[ig] 

steepness  of  the  water  hammer  pressure  pulse,  ETBFCT,  a 
vectorized  continuously  rezonable  flux-corrected  transport^^^^module, 
was  utilized  in  Lagranglai  mode.  Satisfactory  results  were  obtained 
even  on  runs  of  i0,000  to  20,000  timesteps,  provided  r^^  was  taken  to 
be  of  order  5“10  times  the  final  radius  r^.  If  r^  was  larger  than 
this,  however,  grid  point  bunching  occurred,  resulting  in  inordina- 
tely short  timesteps  and  unphysical  profiles. 

The  grid  was  laid  down  initially  so  that  the  spacing  would 
have  become  uniform  if  the  liner  had  imploded  incompressibly  to  the 
origin.  This  prescription  yielded  a closely  gridded  region  near  r^ 
and  large  spacing  near  as  desired.  Grid  point  bunching  was 
alleviated  somewhat  by  filtering  the  flow  velocity  when  it  was  used 
to  advance  the  grid  points  (but  not  in  advancing  the  fluid  equations): 

4.1 

’'j  “ V “j 


where  Uj  was  calculated  using  the  f inlte-dif ferenced  form  of 

^ r I;  ('^2) 


20 


The  best  results  were  obtained  with  D = 0.5(6r)^,  where  6r  is  the 
.variable)  mesh  spacing.  Use  of  instead  of  u.  also  greatly 
reduced  oscillations  in  the  grid  motions.  This  is  particularly 
important  in  the  neighborhood  of  r = r^,  because  the  finite- 
difference  form  of  the  acceleration  of  the  inner  surface  [Eq.  (25 )j 
was  monitored  as  a stability  diagnostic.  At  the  outset  of  the 
calculation  g(t)  displayed  drastic  "ringing"  in  these  grid  points 
for  all  but  the  thickest  liners.  This  was  attributed  to  discrepan- 
cies between  the  initial  profiles  used  and  the  exact  solutions  they 
approximated.  The  ringing  (which  showed  up  in  g but  not  in  r^  or 
r^)  damped  away  quickly,  whereupon  g varied  smoothly  until  after 
turnaround,  when  ringing  again  ensued. 

Aside  from  the  high  frequency  oscillations  superposed  on  g, 
which  were  clearly  of  numerical  origin,  all  the  results  were  sub- 
stantially independent  of  the  number  of  mesh  points,  chosen  to  vary 
between  25  and  200,  and  the  choice  of  r^,  except  as  noted  above. 

The  number  of  mesh  points  was  usually  taken  to  be  the  largest  number 
which  would  permit  running  all  the  cases  with  a given  value  of  M 
in  ~ 10,000  cycles  or  less. 

Conservation  of  mass  and  angular  momentum  were  built  into  the 

s code,  and  were  satisfied  to  within  roundoff,  but  conservation  of 

1; 

■ energy  was  subject  to  truncation  error.  Up  to  turnaround, 

(i  was  constant  to  ^ 2^  in  all  cases.  After  the  water  hammer  was 

I. 

i launched,  however,  W decreased  by  as  much  as  10-15^*  This  was 

ICO  c 

attributed  to  entropy  production,  that  is,  to  the  tendency  of  the 
outgoing  wave  to  steepen  and  form  a shock.  As  noted  above,  our 
equation  of  state  is  isentropic  and  cannot  describe  shocks.  In 
spite  of  this,  we  feel  that  the  overall  wave  motion  is  calculated 
! approximately  correctly,  until  it  reflects  off  the  free  surface  of 

I the  outer  radius. 

IWe  have  dealt  with  the  numerical  techniques  we  used  and  the 
problems  associated  with  them  in  unusually  fine  detail  for  two 
reasons.  First,  the  use  of  an  explicit  code  to  solve  compressible 
fluid  equation  in  the  nearly  Incompressible  limit  is  very  difficult. 


Ik. 


) 

I 

! 

> 

; and  the  answers  obtained  need  to  be  examined  critically.  Second, 

we  have  developed  no  analytical  theory  which  reproduces  the  systema- 
tic content  of  our  results,  and  consequently  cannot  resort  to  theory 
I to  check  or  corroborate  them.  In  spite  of  this,  on  the  basis  of 

considerable  effort  spent  in  code  development  and  numerical  tests, 

r 

j we  feel  confident  in  their  accuracy. 

I III.  Numerical  Results 

t 

A.  Energetics 

[ Using  the  techniques  described  in  the  previous  section  for  the 

I "slow  to  intermediate"  range  of  liner  speeds,  M = 0.01  - 0.2,  we 

obtained  the  results  shown  in  Figs.  2-5*  Curves  are  plotted  for 

I 

1 five  cases,  denoted  by  A through  E.  Cases  A-D  have  a payload  gas 

I with  V = 5/3>  and  M = 0.2,  0.1,  0.05  and  0.01,  respectively,  while  | 

1 for  case  E,  Y = 2 and  M = 0.1.  The  broken  portions  of  the  traces  j 

I represent  interpolations  between  the  calculation  and  the  incompres-  ’ 

I sible  (C^^O)  limit,  in  the  range  where  the  explicit  finite  difference  I 

I methods  used  become  erratic  or  prohibitively  expensive.  The  | 

remainder  of  the  curves  (solid  traces)  are  smooth  fits  through  the 
numerical  results. 

Figure  2 shows  that,  except  for  M « 1 (case  D)  and  « 1, 

the  fraction  of  energy  residing  in  liner  internal  compression  is 

roughly  independent  of  M,  and  y.  A single  curve  fits  almost  all 

the  points  well;  the  flat  part  corresponds  to  ^i^nt^'^tot  **  As 

■*  0,  the  liners  get  thinner,  the  pressure  profiles  across  them  I 

become  flatter,  and  W,  0.  This  requires  the  continuation 

' int  tot  ^ 

(broken  trace)  to  turn  down  as  shown.  For  < 0.05 j the  values  so 

obtained  agree  to  within  a factor  of  two  with  those  calculated  by  | 

Robson^ I 
Case  D evidently  disproves  the  universality  of  this  curve  and 
suggests  the  existence  of  a nontrivial  dependence  on  M.  There  seems 
to  be  little  need  to  pin  this  down,  however,  as  all  points  on  curve  D 
but  the  ones  with  C£  < 0.1  refer  to  extraordinary  thick  liners 
(b/a^  ^ 500;  see  Fig.  8 below). 
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A reduction  in  payload  compression  efficiency  of  order  15^  or 

less  due  to  liner  compressibility  is  quite  acceptable.  However, 

this  is  not  the  whole  story.  The  radial  (translational)  portion  of 

the  kinetic  energy  is  considerably  larger  and  increases  with  both  M 

and  C^.  Together,  W.  and  W , amount  to  about  half  of  ^ when 
f ® ’ rnt  rad  tot 

^ 0.5.  In  neither  plot  does  the  dependence  on  y (represented  by 
the  differences  between  cases  B and  E)  appear  very  significant.  For 
experimentally  realizable  payloads  consisting  of  a mixture  of  plasma 
and  flux,  the  effective  value  of  Y can  be  taken  constant,  approxima- 
tely equal  to  1.3^^^^. 

Figure  shows  that  the  fraction  of  the  total  energy  required 
for  rotational  stabilization  decreases  as  increases.  All  of  the 
traces  have  been  extended  to  intercept  the  ordinate  at  1-Y 
corresponding  to  the  intercepts  at  Y ^ in  Fig.  5«  This  is  in  accord 
with  the  discussion  following  Eq.  (38)  above,  where  it  was  noted 
that  for  an  incompressible  liner,  ^ "*  y-l  in  the  thick-liner 

limit.  Of  course,  “ 0 refers  to  a liner  of  infinitesimal  thick- 
ness, for  which  e=0  and  6 = W /W  =1.  The  crossover  (b~l)  to 
’ rot  tot 

the  thin-liner  case  occurs’ at  such  low  final  pressures,  however,  that 
no  error  is  made  for  finite  in  ignoring  it.  Practically  speaking, 
all  liners  of  interest  are  thick  at  turnaround. 

Because  of  the  Y-dependence  in  the  incompressible  limit,  cases 
B and  E are  necessarily  different  at  low  It  is  seen  from  Fig. 

that  they  also  differ  for  larger  C^.  The  dependence  on  M is  also 
quite  noticeable.  Evidently  there  is  an  energetic  penalty  for  going 
to  fast  implosions  ascribable  to  both  the  radial  and  rotational 
components  of  the  kinetic  energy. 

In  all  cases,  marginal  stability  occurred  before  turnaround. 

The  net  acceleration  g(t)  [Eq.  (25)]  always  climbed  to  a maximum, 
then  dropped  away  as  compressibility  effects  flattened  the  shape  of 
the  trajectory  r(t),  reducing  d^r/dt^.  Because  the  difference  be- 
tween the  inertial  and  centrifugal  terms  in  g is  not  a strong  function 
of  radius  near  turnaround,  the  relation  between  C^,  b and  (the 

amount  of  rotation  needed  for  marginal  stabilization)  is  almost  the 
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same  as  that  for  incompressible  liners,  Eq,  54.  The  discrepancy 
was  ^ for  all  cases  with  b ^ 100,  that  is,  roughly  speaking,  for 
M > 0.01  and  Cf  ^ 0.5. 

figure  5 shows  the  compression  efficiency  e = ^pay^^tot  ^ 
function  of  M and  The  curves  drop  off  with  increasing  very 

sharply  at  first,  then  more  smoothly.  Efficiency  improves  with 
increasing  M,  varying  by  a factor  of  ~ 3 for  a given  value  of 
over  the  range  of  Mach  numbers  investigated.  The  efficiency  is 
significantly  greater  for  Y = 5/5  then  for  y = 2,  the  difference 
being  approximately  twice  the  incompressible  value  of  20^  for  most  of 
the  range  of 

If  we  plot  e ^ against  M and  the  results  can  be  closely 
fitted  by  straight  lines  through  the  incompressible  intercept.  Thus 
we  have  as  a good  approximation 

. yCi  + h(M/Y)C£],  (75) 

where  h is  evidently  an  decreasing  function  of  M and  an  Increasing 
function  of  y [cf.Eq.  (11)  in  the  Introduction,  the  fit  obtained  by 
Schaffer  for  a non-rotating  liner. 1 We  have  no  adequate  model 
to  explain  this  simple  dependence  on  but  the  following  heuristic 
argument  may  at  least  make  it  plausible. 

As  we  have  seen,  treating  the  effects  of  compressibility 
changes  and  only  slightly  from  the  values  found  for 

incompressible  liners,  and  the  changes  have  opposite  signs.  Hence 
Wpay  + roughly  the  same  in  both  cases,  and  we  can  modify 

the  thick-liner  limit  of  Eq.  (58)  to  read 


or 


W + W . 
pay  rad 


W 


tot 


1 

Y ^ 


W 


tot 


W 


pay 


Y(W  +W  ,) 
pay  rad 

W 

pay 


(7^) 

(75) 
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As  we  noted  in  the  Introduction,  the  finite  speed  of  sound  makes  the 
liner  act  as  though  only  a portion  with  outer  radius  r*  < r^  is 
effective  in  compressing  the  payload.  Thus  Eq.  (75)  can  be  written 


.-1 


where 


I 


r,  + K ct. 
£ 1 


Here  t is  the  dwell  time  of  the  liner, 
t - K^r^/u^^^, 


(76) 


(77) 


(78) 


and  are  factors  of  order  unity,  and  the  effective  velocity 
Ugj^  is  determined  from  expression  (5I)  for  the  radial  kinetic 
energy  by 


2 W 


u 


rad 


2 W 


rad 


eff 


TTpr^2  m TTpr^^  im.  {x^lv^) 


(79) 


Hence  by  (75) 


u 


2 _ 


eff 


2 (L-.V) 

eY  n r ^p  ^ (r  ^lx^\ 
1 2 f ' 


« c^ 


2 C. 


Y(Y-l)  e ^ {x^tx^) 


(80) 


Combining  Eqs.  (76-78)  and  (80)  in  the  limit  of  large  where 


r*  - r^  < < r^,  there  results 


,y , \ {x^/x^)  Cf 


2K^K^(Y-1) 


(81) 


While  x^lx^  does  depend  on  the  dependence  in  e"^  is  logarithi- 
mically  weaker  and  can  be  omitted.  Interpolation  between  Eq.  (8I) 
and  the  0 limit  e“^  » Y now  yields  the  desired  form  (75)* 
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Turning  to  the  question  of  thermonuclear  efficiency,  we 
calculate  the  yield  from  Eq.  (70)  according  to 

Y = J dt  (t)  tt  r^2  (t).  (82) 

•00 

[The  actual  limits  of  the  integral  are  from  t » 0 to  approximately 
the  time  when  the  calculation  stops,  i.e.,  when  r is  1 to  2 times 
r (O).]  Using  the  fact  that  is  practically  constant  as 

a function  of  temperature,  we  rewrite  Eq.  (82)  in  the  form 

Y - c^  f cdt/r^  (nT/pc^  )2(  p^^/T2)n  (85) 

The  Q = given  by 

Q - pc  r^  (2/b^^tf)  S cdt/r^  (nT/pc2)2(  P^^^/T^)  (r^/nr^)^,  > 3^) 

where  the  right  hand  side  is  a scaling  factor  Pcr^  multiplying 
an  expression  which  depends  only  on  dimensionless  trajectory 
parameters. 

In  Fig.  (7),  Q/pcr^  is  plotted  against  for  the  cases  A-D 
defined  previously.  It  is  significant  that,  although  Q increases 
noticeably  with  M,  there  is  no  dependence  to  speak  of  on  C£* 

This  may  seem  surprising  in  view  of  the  results  of  Fig.  5^  which 
showed  that  relatively  little  energy  gets  into  the  payload  for 
large  The  explanation  lies  in  the  fact  that  the  dwell  period, 

when  r^  ~ r^  and  is  near  its  maximum,  is  extended  due  to  tamping 
of  the  inner  portion  of  the  liner  by  the  outer  portion.  Figure  8 
shows  how  this  makes  r^(t)  much  flatter  during  the  rebound.  The 
resulting  enhancement  in  the  thermonuclear  burn,  which  is  quadratic 
in  T,  is  even  more  pronounced.  The  beneficial  and  adverse  effects  of 
liner  compressibility  appear  in  this  sense  roughly  to  cancel. 

Another  dimensionless  parameter  which  can  be  used  in  place  of  M 
to  characterize  the  liner  trajectory  is  the  radial  compression  ratio 
■ r^(0)/r^.  For  a liquid  liner  driven  by  pistons,  high  efficiency 
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dictates  that  there  should  be  as  little  fluid  as  possible  in  the 
form  of  an  initial  buffer  layer  on  the  chamber  wall  and  in  the  ducts 
after  the  piston  stroke  is  finished.  To  the  extent  this  is 
achieved,  then,  we  can  set 

r^(0)  rg(t^)  « r^,  (85)  ^ 

where  r is  the  radius  of  the  chamber.  Thus 
w 

r^(t^)/r^  ^ b^.  (86) 

1 

Figure  9 displays  computed  values  of  (i.e.,  b^)  as  a function  of 
and  M for  cases  A-D.  Note  that  the  requirement  of  rotational 
stabilization  results  in  very  thick  (b^  > lOO)  liners  for  all 
pressures  except  0.1  when  M * 0.05. 

Using  Fig.  9 to  determine  M as  a function  of  and  enables 
us  to  plot  Q/pcr^  vs,  the  latter  quantities  (Fig.  lO).  For  a fixed 
compression  ratio,  Q vanishes  as  "*  0.  Q increases  rapidly  at 
first,  then  slowly.  For  a given  final  pressure  Q is  a decreasing 
function  of  a^;  at  large  however,  this  dependence  is  weak. 

From  Fig.  10,  it  is  evident  that  the  thermonuclear  gain 
increases  monotonically  with  increasing  relative  payload  pressure 
The  basic  hydrodynamic  results  for  compressible  liquid  liner 
implosions  thus  do  not  by  themselves  indicate  an  optimum  operating 
value  from  At  fixed  compression  ratio,  higher  values  of 

appear  to  lead  to  smaller  values  of  r^  for  desired  thermonuclear 
gains,  i.e.,  to  smaller  system  sizes.  The  total  system  size,  how- 
ever, will  depend  not  only  on  the  outer  radius  of  the  liner,  0.^r^, 
but  also  on  the  thickness  of  the  pressure  vessel  required  to  contain 
the  liquid  liner  during  implosion  and  reexpansion.  In  an  idealized 
way,  this  consideration  of  the  engineering  requirements  of  the  system 
will  increase  the  total  system  radius  by  a factor  (1  + p^/S),  where 
Pjj  is  the  peak  pressure  of  the  gas  driving  the  liner  implosion  and  S 
is  the  allowable  mechanical  stress.  The  peak  pressure  will  depend 
on  the  particular  driving  mechanism  but  can  be  scaled  by  the  work 
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done  at  constant  pressure  to  obtain  the  desired  payload  pressure: 

CfPC^  ^ 

Pd  * (Y-l)a^^  ' e(Cf,a^)  " 

where  K is  a proportionality  constant  to  account  for  nonconstant 
driving  pressure  during  the  implosion.  The  total  radius  r^,  of  the 
system  may  then  be  estimated  as 

r^  - r^j(l+Pjj/S)  = r^^(l+E  Q^/ea.^^)  (88) 

^ \ Pc^ 

where  Z = A— g—  is  a parameter  based  on  the  properties  of  the  liner, 

payload  and  drive  mechanism,  and  X » K/(y-1). 

The  results  of  Fig.  10  can  be  used  to  obtain  r^  for  a desired 

value  of  Q in  terms  of  mechanical  properties  and  the  peak  payload 

pressure  for  fixed  values  of  a-.  In  Fig.  11,  we  plot 

(Xc/S)Q/r^  versus  for  * 50  two  values  of  E * 50  and  600, 

(corresponding  roughly  to  constant  pressure  drive  with  sodium- 

potassium  and  lead  liners,  respectively,  and  a mechanical  stress  of 

about  10®  Pa).  A clear  maximum  is  evident  for  E = 600,  occurring 

at  “ 0.1,  while  for  E =*  50  the  maximum  is  broader  aid  occurs  at  about 

* 0.^5*  Dividing  the  scaled  values  for  Q displayed  in  Fig.  11 

by  e(C-),  scaled  values  for  the  thermonuclear  energy  gain  Q 
^ P 

relative  to  the  peak  plasma  energy  may  be  obtained  as  shown  in 

Fig.  12.  Note  that,  as  expected,  the  optimum  operating  points 

shift  to  higher  values  of  Cj:*  for  a "scientific  feasibility  experi- 

^ * 

ment,  defined  as  Q =1,  with  = 50  and  E = 600,  the  idealized 
total  system  radius  is  found  by  dividing  the  scaled  value  at 
= 0.2  into  Xc/S.  If  c ■ 1.6  X 10®  cm/ sec,  X * 1.^  and  S » 7*3  x 
10®  dyne/cm^  (10,000  psi),  then  the  necessary  total  radius  is 
r^  * 29  cm.  Magnetic  diffusion  at  the  inside  surface  of  the  liner, 
finite  effective  beta  values  for  the  payload,  and  practical 
engineerlTtb  design  considerations  will,  of  course,  increase  the 
system  size.  The  present  calculations  are  merely  to  indicate  the 
manner  in  which  the  optimum  operating  value  of  Cjj  can  be  determined. 
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B.  Water  Hannier 

Figure  (I3)  shows  the  development  of  the  water  hammer  in  a 
rotating  compressible  liquid  liner.  We  observe  that  the  pressure 
profile  displays  the  characteristic  internal  maximum.  Just  after 
turnaround  (Fig.  I3c),  the  payload  pressure  has  begun  to  drop  from 
its  maximum  {Q^  ^ 0.5)  and  a rarefaction  is  developing.  The  rare- 
faction broadens,  while  the  front  of  the  pulse  (the  water  hammer 
proper)  at  r = r steepens  (Figs.  I5d,  e).  These  processes 
continue  until  the  water  hammer  reaches  the  outside  of  the  liner. 

There  the  free  surface  boundary  condition  (p„  * p(r  ) * O)  forces 

^ 2 

the  pulse  to  reflect  as  a rarefaction,  and  the  calculation  is  no 
longer  physically  meaningful. 

We  can,  however,  cancel  out  the  spurious  effects  of  the  free- 

surface  boundary  condition  sufficiently  to  estimate  the  radial 

sidewall  pressure  produced  when  the  pulse  reaches  a solid  boundary. 

We  do  this  by  noting  that  as  the  pulse  propagates  toward  r , the 

2 

value  of  the  pressure  approaches  a 1/r  dependence  increasingly 
closely.  Thus,  writing 


w 


P r 
max  max 


(89) 


and  plotting  p for  successive  instants  until  r is  a few  mesh 
° max 

points  away  from  r , we  find  a number  which  does  not  depend  on  the 
2 

conditions  near  r = r^.  The  dimensionless  form  of  p^,  namely 

is  plotted  in  Fig.  1^4^  against  and  M for  cases  A-C. 
We  observe  that  although  the  cuirves  increase  monotonically  with 
they  become  quite  flat  for  large  values,  especially  in  case  C.  In 
general,  C^/Cf  is  a sizeable  fraction  (typically  ~‘5)  for  small 
Serious  attention  must  be  given  to  design  of  the  experimental 
apparatus  in  order  that  it  be  able  to  tolerate  such  loads. 

Of  even  greater  importance  is  the  force  on  the  endwalls,  since 
they  must  be  built  so  as  to  permit  access  to  the  chamber  and  there- 
fore contain  holes  or  other  structures.  Figure  15  displays  ©(C^^M)^ 
the  dimensionless  peak  endwall  thrust,  calculated  according  to 
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0pc2  TT  r^2  ^ 2n  J p(r)  rdr  (90 ) 


r 

1 

for  the  same  cases  A-C.  Maximum  0 naturally  always  occurs  just  as 
the  water  hammer  reaches  the  outside  of  the  liner;  as  may  be  seen 
from  Fig.  15,  it  increases  strongly  with  for  large,  unlike 

the  sidewall  pressure  Furthermore,  0.is  largest  for  small  Mach 

numbers,  reflecting  the  scaling  by  The  peak  endwall  pressure 

loading  per  unit  outside  circumference,  0/ct.j,  to  a good  approxima- 
tion is  independent  of  M. 


IV.  Discussion 

We  have  calculated  three  related  but  distinct  effects  of  finite 
r compressibility  in  rotationally  stabilized,  cylindrical  liquid 

liners.  These  are  the  reduction  in  compression  ratio,  enhancement 
of  the  dwell  stage  (and  correspondingly,  of  the  thermonuclear  yield), 
and  the  water-hammer-induced  pressure  loading  on  endwalls  and  side- 
walls.  The  calculations  span  the  region  of  M-C^  parameter  space  | 

describing  conditions  which  must  be  attained  in  a fusion  reactor  and  \ 

which  are  achievable  by  imploding  thick  liners.  The  results  are  ; 

' j 

I embodied  in  the  plots  presented  in  Section  III,  especially  Figs.  5~1 } 

\ 9-12  and  lk-15.  I 

[ ; 

I These  results  are  correct  in  detail  only  for  a system  satisfying 

our  model  assumptions,  namely,  cylindrical  symmetry,  no  dissipation  i 

i 

mechanisms,  and  a payload  of  pure  plasma.  Changing  these  assumptions 

can  affect  the  results  in  two  ways.  First,  as  we  have  seen  in  the  I 

] 

discussion  of  the  allowable  total  stress,  the  introduction  of  more 

realistic  assumptions  can  select  an  operating  point  for  the  system,  or  ! 

at  least  define  a region  in  parameter  space  which  is  most  favorable 
from  the  standpoint  of  the  new  considerations.  Second,  new  physical 

assumptions  can  change  the  curves  themselves.  For  example,  inclusion  j 

of  radiation,  end  losses,  thermal  conduction,  etc.,  in  the  plasma  model  ■ 

reduces  both  the  plasma  compression  efficiency  and  the  value  of  Q.  j 

i 
I 
1 
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To  the  extent  that  such  losses  are  small,  we  can  represent 
them  by  multiplicative  efficiency  factors.  Thus  the  effects  of 
magnetic  diffusion  into  the  liner,  incursion  of  liner  vapor 
impurities  into  the  plasma,  and  the  existence  of  a nearly  plasma- 
free  region  of  magnetic  flux  within  the  payload  volume  can  all  be 
represented  in  a single  geometric  factor  We  say  that  the  payload 

energy  W is  distributed  roughly  uniformly  through  a region  of 
radius  r^^  ^ * r^  + A^,  where  A^  is  the  magnetic  skin  depth.  Within 
this  region  the  plasma  occupies  a volume  of  effective  radius  » 

r^  - A^,  where  A^  is  a layer  of  impurities  and  buffer  field.  Then 
the  energy  in  the  plasma  is 


» W [1-2(A  + A )/r,] 

pay  ' 1 2 f 


e,  W 
f pay 


(91) 


Simularly,  the  driving  energy  from  the  source,  is  converted  by 

the  pistons  into  the  total  liner  energy  W^^^  with  an  efficiency 


m 


W..  ^ / W, 
tot  dr 


(92) 


Hence  the  actual  thermonuclear  efficiency  with  inclusion  of  these 
loss  mechanisms  is 


Y/W 


dr 


! Vc't> 


e e e,Q 
m c f p 


e e-  Q 
m f ^ 


(93) 


where  and  Q are  the  compression  efficiency  (Fig.  5)  and  thermo- 
nuclear efficiency  (Fig.  6)  computed  in  this  paper. 

The  present  work  represents  a completed  stage  of  a continuing 
theoretical  program  in  support  of  the  Linus  concept.  It  will  be 
followed  by  calculations  to  be  presented  in  a later  paper,  in  which 
the  inefficiencies  referred  to  in  the  definition  of  — specifi- 
cally, magnetic  diffusion  and  finite  plasma  beta  — are  treated 
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self-conslstently  along  with  compressibility.  Our  eventual 
objective  is  the  realization  of  the  inherently  favorable  character- 
istics afforded  by  controlled  liner  implosion  in  the  design  and  demon- 
stration of  a viable  reactor  configuration. 


1 


52 


Acknowledgements 


We  wish  to  thank  Mrs.  Peggy  Thume  for  her  efforts  in  preparing 
this,  manuscript  for  publication,  and  Dr.  M.  J.  Schaffer  for  dis- 
cussing the  results  of  his  calculations  with  us  in  advance  of 
publication.  Discussion  with  Dr.  A.  E.  Robson  of  the  results  of 
our  calculations  is  also  acknowledged.  This  work  was  supported  by 
the  Office  of  Naval  Research  and  The  Energy  Research  and  Development 
Administration. 


33 


References 


1.  D.  L.  Book,  A,  L.  Cooper,  R.  Ford,  D.  A.  Hammer,  D.  J.  Jenkins, 
A.  E.  Robson  and  P.  J.  Turchi,  Plasma  Physics  and  Controlled 
Nuclear  Fusion  Research  1976.  Vol.  Ill,  p.  507  (Proc.  VI  Inter- 
national Conference  on  Plasma  Physics  and  Controlled  Nuclear 
Fusion  Research,  Berchtesgaden,  FGR,  6-15  October  1976,  IAEA, 
Vienna,  1977)- 

2.  C.  M.  Fowler,  W.  B.  Garn  and  R.  S.  Caird,  J.  Appl.  Phys 

588  (i960). 

3.  H.  Knoepfel  and  F.  Herlach,  Eds.,  Megagauss  Magnetic  Field 
Generation  by  Explosives  and  Related  Experiments.  EUR  2750. e 
(Euratom,  Brussels,  196b). 

4.  A.  D.  Sakharov,  R.  Z.  Lyudaev,  E.  N.  Smirnov,  Yu.  I.  Plyuschev, 
A.  J.  Pavlovskii,  V.  K.  Chernyshev,  E.  A.  Feoktistova,  E.  I. 
Zharinov  and  Yu.  A.  Zysin,  Dokl.  Akad.  Nauk  SSSR  I65,  65  (I965) 
[Sov.  Phys.  Dokl  IQ,  10^5  (1966)]. 

5.  S.  G.  Alikhanov,  V.  G.  Belan,  G.  I.  Budker,  A.  I.  Ivanchenko 
and  0.  N.  Kichigin,  At.  Energia  53^  (1967)  [Sov.  J.  At. 
Energy  23,  1307  (1967)]. 

6.  D.  L.  Book  and  N.  K.  Winsor,  Phys.  Fluids  17^  662  (197^); 

A.  Barcilon,  D.  L.  Book  and  A.  L.  Cooper,  Phys.  Fluids  17,  1707 

(197i^). 

7.  P.  J.  Turchi,  A.  L.  Cooper,  R.  Ford,  and  D.  J.  Jenkins,  Phys. 
Rev.  Lett.  38,  15^  (1976). 

8.  H.  Knoepfel,  Pulsed  High  Magnetic  Fields  (American  Elsevier, 

New  York,  1970),  Chapt.  10. 

9.  J.  P.  Somon,  "L* Equation  d'Etat  des  Solides,"  LGI  Rept  6k/'^ 

{196k). 

lO.  G.  Lehner,  J.  G.  Linhart  and  J.  P.  Somon,  Nucl.  Fus.  ^ 382 


(1984). 

11.  J.  P.  Somon  and  C.  Jablon,  Nucl.  Fus.  I6,  1040  (1978). 

12.  A.  E.  Robson,  NRL  Memorandum  Report  No.  5^72  (1977). 


15- 


15. 


t 

l’ 

i 


i6. 


r 


17. 

18. 

19. 

20. 

21. 

22. 


M.  J.  Schaffer  (to  be  published). 

P.  J.  Turchi,  Proc.  Seventh  Symposium  on  Engineering  Problems 
of  Fusion  Research,  Knoxville,  Tenn. , 25"28  Oct.  1977. 

J.  P.  Somon,  J.  Fluid  Mech.  ^ 769  (I969);  Ref.  3,  p.  67. 

A.  Barcilon,  D.  L.  Book,  J.  P.  Boris,  A.  L.  Cooper,  K.  Hain, 

P.  C.  Liewer,  A.  E.  Robson,  R.  A.  Shanny,  P.  J.  Turchi  and 

N.  K.  Winsor,  Plasma  Physics  and  Controlled  Nuclear  Fusion 
Research  197^.  Vol.  II,  p.  587  (Proc.  V International  Confer- 
ence on  Plasma  Physics  and  Controlled  Nuclear  Fusion  Research, 
Tokyo,  II-I5  November  197^^  IAEA,  Vienna,  1975 J . 

Samuel  Glasstone  and  R.  H.  Lovberg,  Controlled  Thermonuclear 
Reactions  (Krieger  Publ.  Co.,  Huntington,  N.  Y.,  1975)^ 
pp.  20-25. 

R.  Courant,  K.  0.  Friedrichs,  and  H.  Lewy,  Mathematische 
Annalen  100,  32  (1928)  [transl.  IBM  Journal,  March,  I967, 

pp.  215-25^.] 

J.  P.  Boris,  NRL  Memorandum  Report  No.  5257  (197^)* 

J.  P.  Boris  and  D.  L.  Book,  in  Methods  in  Computational 
Physics . Vol.  I6,  ed.  by  John  Killeen,  [Academic  Press,  New 
York,  1976],  p.  85. 

D.  L.  Book,  D.  A.  Hammer  and  P.  J.  Turchi,  Nucl.  Fus.  1^ 

000  (1978). 

J.  W.  Shearer  and  W.  C.  Condi t,  in  Energy  Storage,  Compression 
and  Switching.  Ed.  W.  H.  Bostick,  V.  Nardi  and  0.  S.  F.  Zucker 
(Plenum  Press,  New  York,  1976),  p.  109. 


55 


J 


Fig.  1 — Schematic  of  model  used.  At  some  time  t = 0 the  liner  has 
inner  and  outer  radii  ri(0)  and  r2(0),  respectively  corresponding  to  the 
dashed  circles.  At  turnaround,  these  become  = r^  and  r2(t^)  (sol- 

id circles).  For  a stable  implosion,  pressure  profiles  p(r)  are  qualitatively 
as  shown.  Azimuthal  velocity  v(r)  decreases  as  1/r. 


'^■''o.o  0.5 

Fig.  2 — Ratio  of  internal  (compressional)  liner  energy  to  total  energy 
at  turnaround  af  a function  of  dimensionless  final  pressure  ? f , for  sev- 
eral choices  of  M and  7.  All  cases  are  marginally  stable  against  Ray- 
leigh-Taylor  modes.  As  (no  liner  compression). 


Fig.  3 — Ratio  at  turnaround  of  radial  kinetic  to  total  energy 
vs  for  the  same  cases  as  in  Fig.  2.  Broken  curves  indicate 
extrapolation  of  numerical  results  to  the  incompressible  limit. 
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Fig.  4 — Ratio  at  turnaround  of  rotational  kinetic  to  total  energy 
vs  ff  for  same  cases  as  in  Figs.  2-3.  Note  that  the  7 = 2 incom- 
pressible limit  differs  from  that  of  the  7 = 5/3  cases. 
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Fig.  5 — Ratio  of  payload  (compressed  plasma/flux  configuration) 
to  total  energy  vs.  for  same  cases  as  in  Figs.  2-4 
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Fig.  6 — Ratio  of  total  energy  to  payload  energy  (inverse  of  quantity 
plotted  in  Fig.  5)  vs  f,  showing  straight-line  fits 


Fig.  7 — Thermonuclear  efficiency  (total  yield  divided  by  total  system 
energy),  scaled  by  pcrf,  vs  f f,  for  cases  A-D 


(e)  (f) 


Fig.  13  — Development  of  water  hammer.  Pressure  profile  is  shown  at  initialization 
(a),  shortly  before  turnaround  (b),  immediately  after  turnaround  (c),  and  in  various 
stages  as  the  steep  pressure  pulse  propagates  outward  (d-f).  Time  and  displacement 
are  scaled  by  c and  r^,. 
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